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Abstract 

Hadronic final states with a hard isolated photon are studied using data taken at 
centre-of-mass energies around the mass of the Z° boson with the OPAL detector at LEP. 
The strong coupling as is extracted by comparing data and QCD predictions for event 
shape observables at average reduced centre-of-mass energies ranging from 24 GeV to 
78 GeV, and the energy dependence of is studied. Our results are consistent with 
the running of as predicted by QCD and show that within the uncertainties of our 
analysis event shapes in hadronic Z" decays with hard and isolated photon radiation can be 
described by QCD at reduced centre-of-mass energies. Combining all values from different 
event shape observables and energies gives as{Mz) = 0.1182±0.0015(stat.)±0.0101(syst.). 



(to be submitted to Eur. Phys. J. C) 



The OPAL Collaboration 



G. Abbiendi^, C. Ainslcy'', P.F. Akcsson^, G. Alexander G. Aiiagnostou^, 
K. J. Anderson^ S. AsaP, D.Axen^^ I.Bailey^^, E. Barberio^'^', T.BarillarP, 
R.J.Barlow^^ R.J.Batley^ P.Bechtle^^ T.Behnke^^, K.W. Bell^^ P.J.Bell^ G. Bella^^ 
A.Bellerive^ G. Benelli^ S.Bethke^i, O.BiebeF, O.Boeriu^ P.Bock^°, 
M.Boutemeur^^, S. Braibant^ R.M. Brown^^ H.J. Burckhart'^, S. Campana^, 
P. Capiluppi^ R.K.Camegie^ A.A. Carter^^ J.R. Carter^ C.Y. Changl^ 
D.G.Charlton\ C. Ciocca^ A. Csilling^s, M.Cuffiani^, S.Dado^o, M. Dallavalle^, A. De 

Roeck^, E.A.De Wolf'', K.Desch^^, B.Dienes^^, J. Dubbert^^, E. DuchovnP, 
G.Duckeck^o, LP. Duerdoth^^ E. Etzion^^ F. Fabbri^, P. Ferrari^ F.Fiedler^o, I. Fleck^ 
M.Ford^^ A.Prey^ P. Gagnoni\ J.W. Gary^ C. Geich-Gimbel^ G. Giacomelli^ 
P. Giacomelli^ M. Giunta^, J. Goldberg^", E. Gross^'^ J. Grunhaus^\ M.Gruwe^, 
A. Gupta^ C.Hajdu^s, M.Hamann^^ G.G.Hanson^ A.HareP, M. Hauschild^ 
CM. Hawkes\ R. Hawkings^ G.Herten^, R.D.Heuer^^, J.C. Hill^ D. Horvath^^'^, 
P. Igo-Kemenes^°, K. Ishii^^, H. Jeremie^'', P. Jovanovic^, T.R. Junk^'*, J. Kanzaki^^'", 
D.Karlen^s, K.Kawagoe^^, T. Kawamoto2^ R.K. Keeler^s, R.G. Kelloggl^ 

B. W. Kennedyi^ S.Kluth^i, T. Kobayashi^^, M.Kobel'^ ^ S. Komamiya^^, T.Kramer^^, 

A. Krasznahorkay Jr.^^'^, P. Krieger^'', J. von Krogh^'^, T. KuhP^, M. Kupper^^, 
G.D. Lafferty^^, H. Landsman^°, D.Lanske^^, D. Lellouch^''^, J. Letts", L. Levinson^^, 
J.Lillich^ S.L.Lloydl^ F.K. Loebingerl^ J.Lu^^.i^ A.Ludwig^'*, J.Ludwig^ W. Mader^'*, 
S.Marcellini^ A.J. Martin^^ T. Mashimo^^ P. Mattig"', J.McKenna^^, 
R.A.McPherson^s, F. Meyers'^, W. Menges^^, F.S. Merritt^, H.Mes^'", N.Meyer^^, 

A.Michelini2, S.Mihara^^, G. Mikenberg^^, D.J.Miller^^ W. Mohr^ T.Mori^^, 
A.Mutter^ K. Nagai^^^ L Nakamura^^.^, H. Nanjo22, H.A.NeaF^ S.W. O'Neale^'*, 
A.Oh^ M.J.Oreglia^ S. Orito^^,*^ CPahl^^, G. Pasztor^ ^ J.R.Paterl^ J.E.Pilcher^ 
J.Pinfold^^, D.E.Plane^, O.Pooth^^^ M. Przybycien^'", A.Quadt^^, K. Rabbertz^'^ 

C. Rembser^, P.RenkeF, J.M. Roney^^ A.M. Rossi^ Y.Rozen^", K.Runge^ K. Sachs^ 

T.SaekP, E.K.G. Sarkisyan^'-'', A.D. Schaile^°, O.Schaile^°, P. Scharff-Hansen^ 
J. Schieck^^, T. Schorner-Sadenius^'^, M.Schroder'', M. Schumacher^, R. Seuster^^''^, 

T.G.Shears^'\ B.C. Shen^'*, P. Sherwoodl^ A. Skujal^ A.M.Smith^ R. Sobie^^, 
S. Soldne^-Remboldl^ F.Spano^-^ A.StahF^ D.Stromi^, R. Strohmer^o, S.Tarem^o, 
M.Tasevsky^''^, R.Teuscher^, M.A. Thomson^ E. Torrence^^, D.Toya^^, L Trigger^''", 
Z. TrocsanyP'^, E. Tsur^^ M.F. Turner- Watson\ LUeda^^, B. UjvarP''^, C.F. Vollmer^o, 
P. Vannerem^, R. Vertesi^^'*^, M. Verzocchi^^, H. Voss'^''?, J. Vossebeld C.P. Ward^, 
D.R.Ward^ P.M.WatkW, A.T.Watson^ N.K.Watson^, P.S. Wells^ T.Wengler^ 
N.Wermes^, G.W. Wilsonl5'^ J.A.Wilson^, G.WolP, T.R.Wyattl^ S. Yamashita^^, 

D.Zer-Zion'', L.Zivkovic^o 

^School of Physics and Astronomy, University of Birmingham, Birmingham B15 2TT, UK 
^Dipartimento di Fisica delL Universita di Bologna and INFN, 1-40126 Bologna, Italy 
^Physikalisches Institut, Universitat Bonn, D-53115 Bonn, Germany 
^Department of Physics, University of California, Riverside CA 92521, USA 

^Cavendish Laboratory, Cambridge CBS OHE, UK 

^Ottawa-Carleton Institute for Physics, Department of Physics, Carleton University, Ot- 
tawa, Ontario KIS 5B6, Canada 



^CERN, European Organisation for Nuclear Research, CH-1211 Geneva 23, Switzerland 
^Enrico Fermi Institute and Department of Physics, University of Chicago, Chicago IL 
60637, USA 

^Fakultat fiir Physik, Albert-Ludwigs-Universitat Freiburg, D- 79104 Freiburg, Germany 
^°Physikalisches Institut, Universitat Heidelberg, D-69120 Heidelberg, Germany 
^^Indiana University, Department of Physics, Bloomington IN 47405, USA 
^^Queen Mary and Westfield College, University of London, London El 4NS, UK 
^^Technische Hochschule Aachen, III Physikalisches Institut, Sommerfeldstrasse 26-28, D- 
52056 Aachen, Germany 

^"^University College London, London WCIE 6BT, UK 

School of Physics and Astronomy, Schuster Laboratory, The University of Manchester 
M13 9PL, UK 

^^Department of Physics, University of Maryland, College Park, MD 20742, USA 
^'''Laboratoire de Physique Nucleaire, Universite de Montreal, Montreal, Quebec H3C 3J7, 

Canada 

^^University of Oregon, Department of Physics, Eugene OR 97403, USA 
^^Rutherford Appleton Laboratory, Chilton, Didcot, Oxfordshire 0X11 OQX, UK 
^"Department of Physics, Technion-Israel Institute of Technology, Haifa 32000, Israel 
^^Department of Physics and Astronomy, Tel Aviv University, Tel Aviv 69978, Israel 
^^International Centre for Elementary Particle Physics and Department of Physics, Uni- 
versity of Tokyo, Tokyo 113-0033, and Kobe University, Kobe 657-8501, Japan 
^^Particle Physics Department, Weizmann Institute of Science, Rehovot 76100, Israel 
^^Universitat Hamburg/DESY, Institut fiir Experimentalphysik, Notkestrasse 85, D-22607 
Hamburg, Germany 

^^University of Victoria, Department of Physics, P O Box 3055, Victoria BC V8W 3P6, 

Canada 

^^University of British Columbia, Department of Physics, Vancouver BC V6T IZl, Canada 
^^University of Alberta, Department of Physics, Edmonton AB T6G 2J1, Canada 
^^Research Institute for Particle and Nuclear Physics, H-1525 Budapest, P O Box 49, 

Hungary 

^^Institutc of Nuclear Research, H-4001 Debrecen, P O Box 51, Hungary 
^"Ludwig-Maximilians-Universitat Miinchen, Sektion Physik, Am Coulombwall 1, D-85748 
Garching, Germany 

''^Max-PIanck-lnstitut fiir Physik, Fohringer Ring 6, D-80805 Miinchen, Germany 
^^Yale University, Department of Physics, New Haven, CT 06520, USA 



" and at TRIUMF, Vancouver, Canada V6T 2A3 

now at University of Alberta 

and Institute of Nuclear Research, Debrecen, Hungary 
^ now at Institute of Physics, Academy of Sciences of the Czech Republic 18221 Prague, 
Czech Republic 

^ and Department of Experimental Physics, University of Debrecen, Hungary 
f and MPI Miinchen 

^ and Research Institute for Particle and Nuclear Physics, Budapest, Hungary 

now at University of Liverpool, Dept of Physics, Liverpool L69 3BX, U.K. 
* now at Dept. Physics, University of lUinois at Urbana-Champaign, U.S.A. 
^ and The University of Manchester, M13 9PL, United Kingdom 



^ now at University of Kansas, Dept of Physics and Astronomy, Lawrence, KS 66045, 
U.S.A. 

' now at University of Toronto, Dept of Physics, Toronto, Canada 

current address Bergische Univcrsitat, Wuppertal, Germany 
" now at University of Mining and Metallurgy, Cracow, Poland 
° now at University of California, San Diego, U.S.A. 
P now at The University of Melbourne, Victoria, Austraha 
^ now at IPHE Universite de Lausanne, CH-1015 Lausanne, Switzerland 

now at lEKP Universitat Karlsruhe, Germany 

* now at University of Antwerpen, Physics Department, B-2610 Antwerpcn, Belgium; sup- 
ported by Interuniversity Attraction Poles Programme - Belgian Science Policy 

* now at Technische Universitat, Dresden, Germany 

" and High Energy Accelerator Research Organisation (KEK), Tsukuba, Ibaraki, Japan 
" now at University of Pennsylvania, Philadelphia, Pennsylvania, USA 

now at TRIUMF, Vancouver, Canada 
^ now at Columbia University 
y now at CERN 
^ now at DESY 

* Deceased 



1 Introduction 



In the theory of strong interactions, Quantum Chromodynamics (QCD) [1-3], the strong 
coupUng constant as is predicted to decrease for high energy or short distance reactions: 
a phenomenon known as asymptotic freedom. The values of Us at different energy scales 
have been measured in e+e^ reactions with different centre-of-mass (cms) energies ranging 
from 35 to 209 GeV and confirm the prediction [4-11]. 

Assuming that photons emitted before or immediately after the TP production do not 
interfere with hard QCD processes, a measurement of at the reduced cms energies, 
V^, of the hadronic system is possible by using radiative multi-hadronic events, i.e. 
e+e~ — >• qq7 events. 

Most photons emitted from the incoming particles before the Z° production (initial 
state radiation, ISR) escape along the beam pipe of the experiment. Measurements of 
cross-sections for hadron production with ISR have been presented by the KLOE and 
BaBar collaborations [12-15]. In e+e^ annihilation to hadrons on the TP peak isolated 
high energy photons observed in the detector are mostly emitted by quarks produced 
in hadronic Z° decays (final state radiation, FSR), because on the Z° peak ISR effects 
are supressed. Measurements of in hadronic events with observed photons have been 
performed by the L3 and DELPHI Collaborations [16,17]. The DELPHI collaboration 
has also measured the mean charged particle multiplicity {nc\i){s') using FSR in [18]. 

When an energetic and isolated photon is emitted in the parton shower the invariant 
mass of the recoiling parton system is taken to set the energy scale for hard QCD processes 
such as gluon radiation. In parton shower models [19-21] the invariant mass of an inter- 
mediate parton or the transverse momentum of a parton branching arc used as ordering 
parameters for the parton shower development. In this picture an energetic and isolated 
photon must be produced at an early stage of the shower evolution and therefore can be 
used to deduce the scale for subsequent QCD processes. The vahdity of this method will 
be studied below using parton shower Monte Carlo programs. 

Here we report on a measurement of as from event shape observables determined from 
the hadronic system in events with observed energetic and isolated photons in the OPAL 
experiment. 

2 Analysis method 

The reduced cms energy, x/s', is defined by 2£'beam-\/l — -E'-y/.E'beam, where is the photon 
energy and £'beam is the beam energy. The flavour mixture of hadronic events in this 
analysis is changed compared to non-radiative Z° decay events. The fraction of up-type 
quarks is larger due to their larger electric charge. However, since the strong interaction is 
blind to quark flavour in the Standard Model, as e.g. demonstrated in [22], the difference 
is not taken into account. 

The determination of ag is based on measurements of event shape observables, which 
are calculated from all particles with momenta Pi in an event: 

Thrust T. The thrust T is deflned by the expression 



The thrust axis ut is the direction n which maximises the expression in parentheses. 

— * 

A plane through the origin and perpendicular to ut divides the event into two 
hemispheres Hi and H2. 

Heavy Jet Mass Mh- The hemisphere invariant masses are calculated using the parti- 
cles in the two hemispheres Hi and H2. We define Mh as the heavier mass, divided 
by y/s . 

Jet Broadening variables Bt and Bwr. These are defined by computing the quantity 

B,= (^B!0^] (2) 



for each of the two event hemispheres, H^, defined above. The two observables are 
defined by 

Bt^Bi + B2 and Bw = max(Si, B2) (3) 
where Bt is the total and Bw is the wide jet broadening. 

C-parameter C. The linear momentum tensor G"'^ is defined by 

Q^ff ^ EiPTpf/IPil ^ ^^^^1^2,3. (4) 

E,- \Pj\ 

The three eigenvalues Xj of this tensor define C with 

C = 3(AiA2 + A2A3 + A3Ai). (5) 



Transition value y^. This observable is given by the value of i/cut in the Durham algo- 
rithm where the number of jets in an event changes from two to three. 

In order to verify that using hadronic Z° decays with hard and isolated final state 
radiation allows one to extract as at a reduced scale y/s' we employ simulated events. We 
use the Monte Carlo simulation programs JETSET version 7.4 [19], HERWIG version 5.9 [20] 
and ARIADNE version 4.08 [21], which have different implementations of the parton shower 
algorithms including simulation of FSR. One sample contains hadronic Z'^ decays with 
FSR and ISR (375 k events) while the other samples are generated at lower cms energies 
without ISR (500 k events each). 

We consider the generated events after the parton shower has stopped (parton-level) 
and calculate event shape observables using the remaining partons. The effective cms 
energy \^ is calculated from the parton four-momenta excluding any final state photons 
and the events arc boosted into the cms system of the partons. The samples are binned 
according to the energy -Efsr of any FSR in intervals of 5 GeV width for EpsR > 10 GeV. 

We observe good agreement between the corresponding distributions obtained from 
the Z° sample with FSR and the lower energy samples. For example. Figure 1 shows 
distributions of the event shape observables 1 — T and Mh for two samples with y/s' — 40 
and 70 GeV. We conclude that within the approximations made in the parton shower 
algorithms, hadronic Z° decays with hard and isolated final state radiation can be used 
to extract measurements of as at reduced scales \/s'. 



3 The OPAL Detector and Event Simulation 



The OPAL detector operated at the LEP e+e" colhder at CERN from 1989 to 2000. 
A detailed description of the detector can be found in [23]. We describe briefly the 
important parts of the detector for this study. In the OPAL coordinate system, the x 
axis was horizontal and pointed approximately towards the centre of LEP, the y axis was 
normal to the z-x plane , and the z axis was in the e~ beam direction. The polar angle, 
9, was measured from the z axis, and the azimuthal angle, 0, from the x axis about the 
z axis. 

The central detector measured the momentum of charged particles and consisted of a 
system of cylindrical drift chambers which lay within an axial magnetic field of 0.435 T. 
The momenta p^y of tracks in the x-y plane were measured with a precision of (Jp/pxy ~ 
0.02% © 0.0015 -p^j/fGeV/c] [24]. 

The electromagnetic calorimeters completely covered the azimuthal range for polar 
angles satisfying | cos^^| < 0.98. The barrel electromagnetic calorimeter covered the polar 
angle range | cos 6*1 < 0.82, and consisted of a barrel of 9440 lead glass blocks oriented 
so that they nearly pointed to the interaction region. The two endcaps were each made 
of 1132 lead glass blocks, aligned along the 2;-axis. Each lead glass block in the barrel 
electromagnetic calorimeter was 10x10 cm^ in cross section, which corresponds to an 
angular region of approximately 40 x 40 mrad^. The intrinsic energy resolution was 
ueIE = 0.2% © 6.3%/^E[GeV] [23]. 

Most electromagnetic showers were initiated before the lead glass mainly because of 
the coil and pressure vessel in front of the calorimeter. An electromagnetic presampler 
made of limited streamer tubes measured the shower position. The barrel presampler 
covered the polar angle range | cos^| < 0.81 and its angular resolution for photons was 
approximately 2 mrad. 

JETSET version 7.4 was used to simulate e+e~ — > qq events, with HERWIG version 5.9 
and ARIADNE version 4.08 used as alternatives. Parameters controUing the hadronisation 
of quarks and gluons were tuned to OPAL LEP 1 data as described in [25,26]. We used 
HERWIG version 5.8 [20], PHOJET version 1.05c [27,28] and VERMASEREN version 1.01 [29] for 
two-photon interactions and KORALZ version 4.02 [30] for e^e^ — > t'^t" events. Generated 
events were processed through a full simulation of the OPAL detector [31] and the same 
event analysis chain was applied to the simulated events as to the data. 4,000,000 fully 
simulated events were generated by JETSET, 200,000 events, 1,000,000 events and 55,000 
events were generated by HERWIG, PHOJET and VERMASEREN while 800,000 events were 
generated by KORALZ. 

4 Event Selection 

4.1 Hadronic Event Selection 

This study is based on a sample of 3 million hadronic Z° decays selected as described 
in [32] from the data accumulated between 1992 and 1995 at cms energy of 91.2 GeV. 
We required that the central detector and the electromagnetic calorimeter were fully 
operational. 

For this study, we apply stringent cuts on tracks and clusters and further cuts on 
hadronic events. The clusters in the electromagnetic calorimeter are required to have 
a minimum energy of 100 MeV in the barrel and 250 MeV in the endcap. Tracks are 



required to have transverse momentum pt > 150 MeV/c with respect to the beam axis, 
at least 40 reconstructed points in the jet chamber, at the point of closest approach a 
distance between the track and the nominal vertex do < 2 cm in the r-0 plane and 
Zq < 25 cm in the z direction. We require at least five such tracks to reduce background 
from e+e" t^t~ and 77 qq events. The polar angle of the thrust axis is required 
to satisfy | cos6't| < 0.9, to ensure that events are well contained in the OPAL detector. 
After these cuts, a data sample of 2.4 x 10^ events remains. 



4.2 Isolated Photon Selection 
4.2.1 Isolation Cuts 

Isolated photons are selected in these hadronic events as follows. Electromagnetic clusters 
with an energy E-^c > 10 GeV are chosen in order to suppress background from soft 
photons coming from the decay of mesons. Accordingly, our signal event is defined as 
an e''"e~ — > qq event with an ISR or FSR photon with energy greater than 10 GeV. We 
use electromagnetic clusters in the polar angle region |cos^ec| < 0-72 corresponding to 
the barrel of the detector, where there is the least material in front of the lead glass, see 
Figure 2 a). Also, the non-pointing geometry of the endcap electromagnetic calorimeter 
complicates the cluster shape fitting explained below. The number of clusters in the data 
which satisfy Eec > 10 GeV and | cos6Iec| < 0.72 is 1,797,532. According to the Monte 
Carlo simulation, 99.3% of these selected clusters come from non-radiative multi-hadronic 
events. 

The candidate clusters are required to be isolated from any jets, and from other clusters 
and tracks: 

• The angle with respect to the axis of any jet, ctiso, is required to be larger than 25°, 
see Figure 2 b). The jets are reconstructed from tracks and electromagnetic clusters, 
excluding the candidate cluster, using the Durham algorithm [33] with ^cut = 0.005. 

• The sum of the momenta Piso of tracks falling on the calorimeter surface inside a 0.2 
radian cone around the photon candidate is required to be smaller than 0.5 GeV/c 
(Figure 2c)). The total energy deposition in the electromagnetic calorimeter within 
a cone of 0.2 radian around the photon candidate, £'iso, is also required to be less 
than 0.5 GeV (Figure 2 d)). 

After the isolation cuts, 11,265 clusters are retained. The fraction of clusters from non- 
radiative multi-hadronic events is reduced to 52.8%. The background from t~^t~ events 
(two-photon events) is 0.6% (0.01%) [34]. 



4.2.2 Likelihood Photon Selection 

Isolated photon candidates are selected by using a likelihood ratio method [35] with four 
input variables. The first two variables are |cos^ec| and aiso, defined above. Two more 
variables, the cluster shape fit variable S and the distance A between the electromagnetic 
calorimeter cluster and the associated presampler cluster, defined as follows, reduce the 
background from clusters arising from the decays of neutral hadrons. 
The cluster shape fit variable, S, is defined by 

1 y^ (-E'meaSjj -^exp,i) ^g>^ 



-^block j "^measji 



where A^biock is the number of lead glass blocks included in the electromagnetic cluster, 
Emeas,i is the measured energy deposit in the ith block, -Eexp,i is the expected energy 
deposit in the ith block, assuming that the energy is deposited by a single photon, and 
c"meas,i is the Uncertainty in the energy measured by ith electromagnetic calorimeter block. 
Ecxp,i is a function of position and energy of the incident photon based on the simulation 
of the OPAL detector with single photons. The value of S is determined by minimizing 
Equation (6) under variation of the position and energy of the cluster. For a cluster to be 
considered further in the likelihood, preselection cuts are applied: we require the number 
of blocks to be at least two and the value of S after the fit to be smaller than 10. The 
quality of the cluster shape fits depends on the assumed resolution (Jmeas/, this will be 
studied as a systematic uncertainty. 

The variable A measures the distance between the electromagnetic calorimeter cluster 
and the associated presampler cluster, A = max(|A0|, |A^|), with A0 and A^ the angular 
separations between the clusters. 

The distributions of 5* and A are shown in Figures 2 e) and f). The Monte Carlo 
distributions in these figures are normalized according to the luminosity obtained from 
small angle Bhabha events. 

A disagreement between data and Monte Carlo is seen for S and ctiso. The level of 
agreement between data and Monte Carlo for the S distribution is studied with photons 
in radiative muon pair events and tt^s produced in tau pair events. It is confirmed that the 
Monte Carlo adequately reproduces the S distributions [34]. The disagreement between 
data and Monte Carlo for distributions of S and aiso stems from the failure of the Monte 
Carlo generators to correctly predict the rate of isolated neutral hadrons, as explained in 
Section 4.4. In this analysis, the rate of isolated neutral hadrons used in the background 
subtraction is estimated from data by methods described in Section 4.4. 

The likelihood calculation is performed with reference histograms made for seven sub- 
samples, chosen according to the cluster energy. The cut on the likelihood value is chosen 
so as to retain 80% of the signal events. The hkelihood distributions for data and Monte 
Carlo arc shown in Figure 3. It can be seen that the likelihood distributions for signal and 
background events are well separated for each region of electromagnetic cluster energy. 
The correlations between variables arc taken into account using the algorithm described 
in [36]. Electromagnetic clusters which pass the likelihood selection are regarded as pho- 
ton candidates. If more than one candidate is found in the same event the one with the 
highest energy is chosen. 

4.3 Final Data Sample 

Hadronic events with hard isolated photon candidates are divided into seven subsamples 

according to the photon energy for further analysis. Table 1 shows the mean values of 
y/s', the number of data events and the number of background events for each subsample. 

4.4 Background Estimation 

According to the Monte Carlo simulation, the contamination from r pair events is between 
0.5 and 1.0%. The impact of this small number of events is further reduced because the 
value of event shape observables for r pair events are concentrated in the lowest bin of 
the distributions, outside the fitting range, so their effect on the ccg fits is negligible. The 
contribution of two photon processes is less than 0.01% in all subsamples and is ignored. 



As mentioned in [37], the JETSET Monte Carlo fails to reproduce the observed rate 
of isolated electromagnetic clusters, both for isolated photons and isolated 7r°'s. Isolated 
neutral hadrons are the dominant source of background for this analysis, and their rate 
has been estimated from data using the following two methods. 

Firstly, with the likelihood ratio method the observed likelihood distributions in the 
data in bins of photon energy were fitted with a linear combination of the Monte Carlo 
distributions for signal and background events which pass the isolation cuts and hkelihood 
preselection requirements. The overall normalisation of the Monte Carlo distribution is 
fixed to the number of data events. The fit uses a binned maximum likelihood method 
with only the fraction of background events as a free parameter. Figure 3 shows the fit 
results. The values of x^/d.o.f. are between 1.2 and 3.4 for 18 degrees of freedom. 

Secondly, with the isolated tracks method the fraction of background from isolated 
neutral hadrons was estimated from the rates of isolated charged hadrons. We select from 
the data tracks which satisfy the same isolation criteria as the photon candidates. The 
composition of these isolated charged hadrons obtained from JETSET is used to infer the 
rates of charged pions, kaons and protons. When isospin symmetry is assumed, the rates 
of neutral pions, neutral kaons and neutrons can be estimated from the rates of charged 
pions, charged kaons and protons, respectively: 

-RttO — o-^Tri Rk° — o-Rk± Rn — Rp , (7) 

z ^ z 

where Rx is the production rate of particle X. According to JETSET tuned with OPAL 
data, the rate of isospin symmetry violation is 10% for pions and 5% for kaons and protons. 
This is assigned as a systematic uncertainty for the isolated tracks method and combined 
with the statistical uncertainty. 

The neutral hadron background fractions estimated by these two methods are shown 
in Table 1. The statistical errors from the number of data and Monte Carlo events from 
fitting the hkelihood distributions are shown. The results from the two methods are within 
at most three standard deviations of these errors, except in the bin 35 — 40 GeV. 

The standard analysis will use the likelihood ratio method. Any differences in the 
resulting values of Og obtained by using the two background estimate methods will be 
treated as a systematic uncertainty. 

5 Measurement of Event Shape Distributions 

In this analysis event shape observables as defined above in section 2 are calculated 
from tracks and electromagnetic clusters excluding the isolated photon candidate. The 
contributions of electromagnetic clusters originating from charged particles are removed 
by the method described in [38]. 

We evaluate the observables in the cms frame of the hadronic system. The Lorentz 
boost is determined from the energy and angle of the photon candidate. When the four- 
momentum of particles in the hadronic system is calculated, electromagnetic clusters are 
treated as photons with zero mass while tracks of charged particles are treated as hadrons 
with the charged pion mass. 

Distributions of the event shape observables (1 — T) and Mh are shown for two cms 
energies in Figure 4. The remaining background is removed by subtracting the scaled 
Monte Carlo predictions for non-radiative hadronic events and for r pair events using the 



background estimates listed in Table 1. The effects of the experimental resolution and 
acceptance are unfolded using Monte Carlo samples with full detector simulation (detector 
correction). The unfolding is performed bin- by-bin with correction factors rf^* = hi/di, 
where hi represents the value in the ith bin of the event shape distribution of stable 
hadrons in the Monte Carlo simulation, where "hadrons" are defined as particles with a 
mean proper lifetime longer than 3 • 10~^°s. di represents the value in the ith bin of the 
event shape distribution calculated with clusters and tracks obtained from Monte Carlo 
samples with detector simulation after the complete event selection has been applied. We 
refer to the distributions after applying these corrections as data corrected to the hadron 
level. 

The distributions of the event shape observables 1 — T and Mh for data corrected to 
the hadron level and corresponding Monte Carlo predictions are shown in Figure 5^. The 
Monte Carlo samples are generated with cms energies set to the mean value of in 
each subsample. In the production of the Monte Carlo samples ISR and FSR is switched 
off and on, respectively. The predictions from the event generators are consistent with 
the data for all x/i' bins. There is similar agreement between data and event generator 
predictions for the other observables. 

6 Measurement of cxs 

The measurement of is performed by fitting perturbative QCD predictions to the event 
shape distributions corrected to the hadron level for (1 — T), Mh [39], Brp-, By/ [40,41], 
C [42,43] and [33,44,45]. The 0{a1) and NLLA calculations are combined with the 
ln(i?) matching scheme. The effects of hadronisation on event shapes must be taken into 
account in order to perform fitting at the hadron level (hadronisation correction). Pre- 
serving the normalisation in the hadronisation correction is not trivial for low x/s' samples 
because of large hadronisation corrections. The hadronisation correction is applied to the 
integrated (cumulative) theoretical calculation to conserve normalisation as in our pre- 
vious analyses [46,47]. The hadron level predictions are obtained from the cumulative 
theoretical calculation multiplied by a correction factor Rf"^ — Hi/ Pi, where Pi (Hi) 
represents the value in the ith bin of the cumulative event shape distribution calculated 
by Monte Carlo simulation without (with) hadronisation. The JETSET Monte Carlo event 
generator is used for our central results, while HERWIG and ARIADNE are considered as 
alternatives for the estimation of systematic uncertainties. 

The fit of the hadron level QCD predictions to the event shape observables uses a 
least method with as{Q) treated as a free parameter. Only statistical uncertainties 
are taken into account in the calculation of x^. When the total number of events is small, 
the differences between the statistical errors counting larger or smaller numbers of events 
than the theoretical prediction can bias the fit result. In order to avoid this bias the value 
of the fitted theoretical distribution is used to calculate the statistical error instead of 
the number of events in each bin of the data distribution. The statistical uncertainty is 
estimated from the fit results derived from 100 Monte Carlo subsamples with the same 
number of events as selected data events. 

The region used in the fit is adjusted such that the background subtraction and the 
detector and hadronisation corrections are small (less than 50%) and uniform in that 

^The values of the six observables at the seven energy points are given in [34] and will be available 
under http : //durpdg . dur . ac . uk/HEPDATA/. 



region. The resulting fit ranges are mainly restricted by the hadronisation corrections. 
The QCD predictions at y/s' — 78 GeV fitted to data after applying the hadronisation 
correction are shown in Figure 6. Good agreement between data and theory is seen. The 
fitted values of ag and their errors for each event shape observable are shown in Tables 2-5. 

6.1 Systematic Uncertainties 

6.1.1 Experimental Uncertainties 

The experimental uncertainty is estimated by adding in quadrature the following contri- 
butions: 

• the difference between the standard result and the result when all clusters and tracks 
are used without correcting for double counting of energy. This variation is sensitive 
to imperfections of the detector simulation. 

• the largest deviation between the standard result and the result when the analysis is 
repeated with tighter selection criteria to eliminate background (standard values in 
brackets): the thrust axis is required to lie in the range | cos^t| < 0.7 (0.9), or the 
cluster shape variable is required to be smaller than 5 (10), or the isolation angle 
from any jet is required to be larger than 35° (25°). 

• the difference between the standard value and the value obtained by repeating the 
analysis with the background fractions estimated from the rate of isolated charged 
hadrons as described in Section 4. 

• the difference between the standard result and the result when the single block 
energy resolution is varied to give the expected the cluster shape fits. This 
check is made, because the values of in the cluster shape fits depend on the 
assumed energy resolution. 

• the maximum difference between the standard result and the result when the fit 
regions are varied. The lower and upper limit of the fitting region are independently 
changed by ±1 bin. 

The tighter selection on | cos^r| and the alternative single block energy resolution of the 
electromagnetic calorimeter yield the largest contributions to the experimental systematic 
uncertainty. The overall resolution and energy scale uncertainty of the electromagnetic 
calorimeter have a neglegible effect on the results of this analysis. 

6.1.2 Hadronisation Uncertainties 

The following variations are performed in order to estimate the hadronisation uncertain- 
ties: 

• the largest of the changes in ag observed when independently varying the hadroni- 
sation parameters h and ctq by ±1 standard deviation about their tuned values in 
JETSET [25]; 

• the change observed when the parton virtuality cut-off parameter Qq is varied by 
±1 standard deviation about its tuned value in JETSET; 



• the change observed when only the hght quarks u, d, s and c are considered at the 
parton level in order to estimate potential quark mass effects; 

• the differences with respect to the standard result when HERWIG or ARIADNE are used 
for the hadronisation correction, rather than JETSET. 

We define the hadronisation correction uncertainty by adding in quadrature the deviation 
when using only light quarks and the larger deviation when using HERWIG or ARIADNE to 
calculate the corrections. These variations are observed to lead to larger differences than 
all other variations, i.e. the main contributions to the hadronisation uncertainties are the 
choice of hadronisation model and the potential effect of quark masses. 

6.1.3 Theoretical Uncertainties 

We fix the renormalisation scale parameter = yu/Q to 1, where fi is the energy scale at 
which the theory is renormalized and Q is the energy scale of the reaction. Although the 
uncertainty on the choice of the value of X/^ gives a large contribution to the systematic 
uncertainty, the means of quantifying this uncertainty is essentially arbitrary. We define 
the scale uncertainty as the larger of the deviations of as when is changed from 1 to 
0.5 or 2.0. 

The 0{al) and NLLA calculations are combined with the ln(R) matching scheme. 
The variation in ^^(v^) due to using different matching schemes is much smaller than 
the renormalisation scale uncertainty [48] , and is not included as an additional theoretical 
systematic uncertainty. 

6.2 Combination of Results 

The values of ag obtained for each observable at each energy are used to study the energy 
dependence of ag and to obtain an overall combined result for as{Mz)- The individual 
values of a.^ as given in Tables 2-5 and shown in Figure 7 are combined taking the cor- 
relations between their statistical and systematic errors into account using the method 
described in [8] . The statistical covariances between results from different observables are 
determined at each energy from 100 Monte Carlo subsamples with the same number of 
events as selected in the data. The experimental systematic uncertainties are assumed 
to be partially correlated, i.e. covij = min((Tj, 0"^)^. The hadronisation and theoretical 
covariances are only added to the diagonal of the total covariance matrix. The correla- 
tions between these uncertainties are considered by repeating the combination procedure 
with different hadronisation corrections (udsc only, HERWIG, ARIADNE) and with different 
renormalisation scale parameters {x^ = 2 and 0.5). The systematic uncertainties for the 
combined value are obtained by repeating the combination for each systematic variation. 
The resulting values of as{\/s') are shown in Table 6 and Figure 8. 

Values of as from individual observables at each energy are combined after evolving 
them to a/s — Mz- In this case the results are statistically uncorrelated. The correlations 
between systematic uncertainties are treated as explained above. The results are given in 
Table 7 and Figure 9. 

We also combine the combined values listed in Table 7 taking into account their 
statistical correlations using the sum of the inverses of the individual statistical covariance 
matrices at each energy point. The result is 



as{Mz) = 0.1182 ± 0.0015(stat.) ± 0.0101(syst.). 



(8) 



and is shown with individual errors in Table 7. Figure 8 shows the evolution of the 
strong coupling using our result. As a crosscheck on the robustness of the combination 
procedure we repeat the combination using the combined results at each energy point 
shown in Tabic 6 or using all individual results and find as{Mz) — 0.1183 ± 0.0103 or 
a^(Mz) = 0.1179 ± 0.0103, respectively. 

Our result is consistent within the statistical and experimental errors with the result 
from OPAL using non-radiative events in LEP 1 data with the same set of observables, 
as{Mz) = 0.1192±0.0002(stat.)±0.0050(syst.) [8]. Our result is also consistent with recent 
combined values [11,49-51] and results from other analyses using radiative events [16,17]. 
Figure 10 compares our result with the LEP 1 value for as from [8] and an average of 
results from L3 using radiative hadronic Z° decays [16]^. 

The combinations of individual observables at different cms energies yield x^/d.o.f. ^ 
1/6. The small values of x^/d.o.f. are due to the conservative treatment of hadronisation 
and theoretical uncertainties. The values of x^/d.o.f. indicate consistency of the individual 
results with the model of the combination including evolution of results at different cms 
energies to Mz before the combination. 

7 Summary 

The strong couphng ag has been measured at reduced cms energies, y/s', ranging from 
20 GeV to 80 GeV using event shape observables derived from the hadronic system in 
radiative hadronic events. 

Fits of 0{af) and NLLA QCD predictions to the six event shape observables 1 — T, 
Mh, Bt, BiY, C and |/£ are performed and values of are obtained for seven values of 
sfs'. Our results are consistent with the running of ag as predicted by QCD. The values 
at each y/s' are evolved to jj, — Mz and combined for each event shape observable. The 
combined value from all event shape observables and values is as{Mz) — 0.1182 ± 
0.0015(stat.) ±0.0101(syst.). 

This result agrees with previous OPAL analyses with non-radiative LEP 1 data, with 
a similar measurement by L3, and with recent world average values, see figure 10. Within 
errors, QCD is consistent with our data sample of events with isolated FSR. Our result 
supports the assumption that the effects of high energy and large angle FSR on event 
shapes in hadronic Z° decays can be effectively described by QCD with a lower effective 
cms energy \fs' . 
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Table 1: The number of selected events and the mean value of y/s' for each y/s' subsample. 
The neutral hadron background fractions estimated by the two methods described in 
Section 4.4 are listed in the columns "Non-rad. MH". 
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Table 2: Values of as and their errors for subsamples £'.^=10-15 GeV (upper) and 15- 
20 GeV (lower). 
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± 0.0114 
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-0.0079 


-0.0063 


-0.0078 


-0.0042 


-0.0072 


-0.0023 


a;^ = 2.0 


0.0101 
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0.0098 


0.0058 
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Total error 


+ 0.0175 
- 0.0164 


+ 0.0186 
- 0.0176 


+ 0.0172 
- 0.0162 


+ 0.0130 
- 0.0124 


+ 0.0183 
- 0.0175 


+ 0.0189 
- 0.0180 



Table 3: Values of as and their errors for subsamples — 20 — 25 GeV (upper) and 
25 -30 GeV (lower). 
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Table 4: Values of as and their errors for subsamples — 30 — 35 GeV (upper) and 
35 -40 GeV (lower). 
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Table 5: Values of ag and their errors for subsample — AO — 45 GeV. 



Vs' [GeV] 


78.1 71.8 65.1 57.6 49.0 38.5 24.4 




0.1153 0.1242 0.1201 0.1296 0.1353 0.1438 0.1496 


Statistical 
Experimental 
Hadronisation 
Theory 


0.0026 0.0037 0.0039 0.0047 0.0053 0.0064 0.0071 
0.0068 0.0036 0.0040 0.0069 0.0039 0.0063 0.0077 
0.0062 0.0065 0.0072 0.0085 0.0100 0.0122 0.0166 
0.0053 0.0067 0.0063 0.0076 0.0086 0.0099 0.0117 



Table 6: Combined values of as{\/s') and their errors from all event shape variables. 





(1-T) Mh Bt Bw C y^^ 


Combined 


asiMz) 


0.1230 0.1187 0.1214 0.1117 0.1195 0.1261 


0.1182 


Statistical 
Experimental 
Hadronisation 
Theory 


0.0028 0.0024 0.0021 0.0021 0.0023 0.0031 
0.0050 0.0054 0.0037 0.0033 0.0040 0.0049 
0.0071 0.0052 0.0080 0.0061 0.0092 0.0105 
0.0076 0.0059 0.0081 0.0049 0.0076 0.0045 


0.0015 
0.0038 
0.0070 
0.0062 



Table 7: Combined values of as{Mi) and their errors from all photon energy subsamples 
for a given observable. The final combined value of as{Mz) is also shown. 






Figure 1: The distributions of event shape observables 1 — T and Mh for non-radiative 
events and radiative hadronic events from the Monte Carfo generators JETSET, HERWIG 
and ARIADNE as indicated befow the figures. The triangles and points show distributions 
obtained from the Z° samples with FSR while the histograms show distributions from 
samples generated at lower energies as shown on the figure. The open triangles and solid 
histogram (solid points and dashed histogram) in each figure correspond to y/s' = 40 
(70) GeV. 



OPAL 



t/3 



B fir 
■3 10 



□ MH 

[1 Radiative MH 
• Data 



<— I 



a) 





100 150 

^un /degrees 





Figure 2: Distributions of each variable used in the isolated photon selection. The error 
bars show the statistical errors. Monte Carlo distributions are normalized to the integrated 
luminosity of the data and the cross section of the process. Arrows in the figures show the 
selected region. Distributions for radiative multi-hadronic events, which are signal events 
in this analysis, are overlaid on distributions for all multi-hadronic events and rr events. 
The distribution of each variable is obtained with the cuts on the preceeding variables 
applied. 
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Figure 3: Photon likelihood distributions. The error bar shows statistical error. The 
Monte Carlo distributions are normalised to the total number of candidates in the data, 
and the neutral hadron background fractions are obtained from the fits described in 
Section 4.4. The arrows indicate the selected regions. 
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Figure 4: Event shape distributions before background subtraction and detector correc- 
tion. Two of the six event shape observables, 1 — T and Mh, are shown for the low 
(38.5 GeV) and high (78.1 GeV) Vs' samples. The histograms show Monte Carlo distri- 
butions. The error bars show the statistical errors. 
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Figure 5: Event shape distributions at the hadron level. The error bars correspond to 
the statistical and experimental uncertainties described in Section 6.1.1. Two of the 
six event shape observables, 1 — T and M/^, are shown for the low (38.5 GeV) and high 
(78.1 GeV) \fs' samples. The small lines on the error bars show the extent of the statistical 
uncertainty. The data points are placed at the centres of the corresponding bins. The 
predictions of JETSET, HERWIG and ARIADNE at the corresponding values are also 
shown as lines. 
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Figure 6: Event shape distributions for data at (v s') = 78.1 GeV and the fitted theoretical 
predictions. The error bars show the statistical errors. The solid lines in the theoretical 
predictions show the regions used in the fit. Three corrections are plotted as "i?corr": 
the detector correction, rf^* (dashed line), the hadronisation correction, Rf"''^ (solid line), 
and the ratio of distributions after and before background subtraction (dotted line). The 
hadronisation correction is shown by the ratio of differential distributions in these figures 
(see text for details). 
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Figure 7: Energy dependence of as for all ^fs' subsamples. The inner error bars show the 
statistical and the outer error bars the total uncertainties. The curves and shaded bands 
show the QCD prediction for the running of as obtained with the corresponding values 
of as{Mi) with total errors from Table 7. 




Figure 8: Combined values of from all event shape observables as shown in Table 6. 
The curve and shaded band show the QCD prediction for the running of as using the 
combined value of as{Mz) with total errors. 
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Figure 9: The values of as{Mz) obtained by combining all y/s' samples as shown in 
Table 7. The inner error bars are statistical, the outer error bars correspond to the total 
uncertainty. The dashed vertical lines and shaded bands show the LEP 1 results from 
OPAL [8] using non-radiative events. 
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Figure 10: Combined values of as{Mz) for all event shape observables and \/^ samples. 
The error bars show total uncertainties. The results from this analyses with radiative 
events, from non-radiative events with OPAL LEP 1 data [8] and from L3 radiative 
events [16] are shown. The PDG [49] value of as{Mz) is also shown as the vertical line, 
with the total uncertainty corresponding to the shaded band. 
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1 Introduction 



In the theory of strong interactions, Quantum Chromodynamics (QCD) [1-3], the strong 
coupUng constant as is predicted to decrease for high energy or short distance reactions: 
a phenomenon known as asymptotic freedom. Values of Ug at different energy scales have 
been measured at PETRA and LEP in e^e" reactions with different centre-of-mass (cms) 
energies ranging from 35 to 209 GeV and confirm the prediction [4-11]. 

Assuming that photons emitted before or immediately after the TP production do not 
interfere with hard QCD processes, a measurement of at the reduced cms energies, 
V^, of the hadronic system is possible by using radiative multi-hadronic events, i.e. 
e+e~ — >• qq7 events. 

Most photons emitted from the incoming particles before the Z° production (initial 
state radiation, ISR) escape along the beam pipe of the experiment. Measurements of 
cross-sections for hadron production with ISR have been presented by the KLOE and 
BaBar collaborations [12-15]. In e+e^ annihilation to hadrons on the TP peak isolated 
high energy photons observed in the detector are mostly emitted by quarks produced in 
hadronic Z° decays (final state radiation, FSR), because on the Z° peak ISR effects are 
suppressed. Measurements of ag in hadronic events with observed photons have been 
performed by the L3 and DELPHI Collaborations [16,17]. The DELPHI collaboration 
has also measured the mean charged particle multiplicity {nc\i){s') using FSR in [18]. 

When an energetic and isolated photon is emitted in the parton shower the invariant 
mass of the recoiling parton system is taken to set the energy scale for hard QCD processes 
such as gluon radiation. In parton shower models [19-21] the invariant mass of an inter- 
mediate parton or the transverse momentum of a parton branching arc used as ordering 
parameters for the parton shower development. In this picture an energetic and isolated 
photon must be produced at an early stage of the shower evolution and therefore can be 
used to deduce the scale for subsequent QCD processes. The vahdity of this method will 
be studied below using parton shower Monte Carlo programs. 

Here we report on a measurement of ag from event shape observables determined from 
the hadronic system in events with observed energetic and isolated photons in the OPAL 
experiment. 

2 Analysis method 

The reduced cms energy, x/s', is defined by 2£'beam-\/l — -E'-y/E'beam, where is the photon 
energy and £'beam is the beam energy. The flavour mixture of hadronic events in this 
analysis is changed compared to non-radiative Z° decay events. The fraction of up-type 
quarks is larger due to their larger electric charge. However, since the strong interaction is 
blind to quark flavour in the Standard Model, as e.g. demonstrated in [22], the difference 
is not taken into account. The effects of massive b quarks on hadronisation corrections 
are considered below as a systematic uncertainty. 

The determination of ag is based on measurements of event shape observables, which 
are calculated from all particles with momenta Pi in an event: 

Thrust T. The thrust T is deflned by the expression 

rj. I Et \Pi ■ n\ \ 



(1) 



The thrust axis ut is the direction n which maximises the expression in parentheses. 

— * 

A plane through the origin and perpendicular to ut divides the event into two 
hemispheres Hi and H2. 

Heavy Jet Mass Mh- The hemisphere invariant masses are calculated using the parti- 
cles in the two hemispheres Hi and H2. We define Mh as the heavier mass, divided 
by y/s . 

Jet Broadening variables Bt and Bwr. These are defined by computing the quantity 

B,= (^B!0^] (2) 



for each of the two event hemispheres, H^, defined above. The two observables are 
defined by 

Bt^Bi + B2 and Bw = max(Si, B2) (3) 
where Bt is the total and Bw is the wide jet broadening. 

C-parameter C. The linear momentum tensor G"'^ is defined by 

Q^ff ^ EiPTpf/IPil ^ ^^^^1^2,3. (4) 

E,- \Pj\ 

The three eigenvalues Xj of this tensor define C with 

C = 3(AiA2 + A2A3 + A3Ai). (5) 



Transition value y^. This observable is given by the value of i/cut in the Durham algo- 
rithm where the number of jets in an event changes from two to three. 

In order to verify that using hadronic Z° decays with hard and isolated final state 
radiation allows one to extract as at a reduced scale y/s' we employ simulated events. We 
use the Monte Carlo simulation programs JETSET version 7.4 [19], HERWIG version 5.9 [20] 
and ARIADNE version 4.08 [21], which have different implementations of the parton shower 
algorithms including simulation of FSR. One sample contains hadronic Z'^ decays with 
FSR and ISR (375 k events) while the other samples are generated at lower cms energies 
without ISR (500 k events each). 

We consider the generated events after the parton shower has stopped (parton-level) 
and calculate event shape observables using the remaining partons. The effective cms 
energy \^ is calculated from the parton four-momenta excluding any final state photons 
and the events arc boosted into the cms system of the partons. The samples are binned 
according to the energy -Efsr of any FSR in intervals of 5 GeV width for EpsR > 10 GeV. 

We observe good agreement between the corresponding distributions obtained from 
the Z° sample with FSR and the lower energy samples. For example. Figure 1 shows 
distributions of the event shape observables 1 — T and Mh for two samples with y/s' — 40 
and 70 GeV. We conclude that within the approximations made in the parton shower 
algorithms, hadronic Z° decays with hard and isolated final state radiation can be used 
to extract measurements of as at reduced scales \/s'. 



3 The OPAL Detector and Event Simulation 



The OPAL detector operated at the LEP e+e" colhder at CERN from 1989 to 2000. 
A detailed description of the detector can be found in [23]. We describe briefly the 
important parts of the detector for this study. In the OPAL coordinate system, the x 
axis was horizontal and pointed approximately towards the centre of LEP, the y axis was 
normal to the z-x plane , and the z axis was in the e~ beam direction. The polar angle, 
9, was measured from the z axis, and the azimuthal angle, 0, from the x axis about the 
z axis. 

The central detector measured the momentum of charged particles and consisted of a 
system of cylindrical drift chambers which lay within an axial magnetic field of 0.435 T. 
The momenta p^y of tracks in the x-y plane were measured with a precision of (Jp/pxy ~ 
0.02% © 0.0015 -p^j/fGeV/c] [24]. 

The electromagnetic calorimeters completely covered the azimuthal range for polar 
angles satisfying | cos^^| < 0.98. The barrel electromagnetic calorimeter covered the polar 
angle range | cos 6*1 < 0.82, and consisted of a barrel of 9440 lead glass blocks oriented 
so that they nearly pointed to the interaction region. The two endcaps were each made 
of 1132 lead glass blocks, aligned along the 2;-axis. Each lead glass block in the barrel 
electromagnetic calorimeter was 10x10 cm^ in cross section, which corresponds to an 
angular region of approximately 40 x 40 mrad^. The intrinsic energy resolution was 
ueIE = 0.2% © 6.3%/^E[GeV] [23]. 

Most electromagnetic showers were initiated before the lead glass mainly because of 
the coil and pressure vessel in front of the calorimeter. An electromagnetic presampler 
made of limited streamer tubes measured the shower position. The barrel presampler 
covered the polar angle range | cos^| < 0.81 and its angular resolution for photons was 
approximately 2 mrad. 

JETSET version 7.4 was used to simulate e+e~ — > qq events, with HERWIG version 5.9 
and ARIADNE version 4.08 used as alternatives. Parameters controUing the hadronisation 
of quarks and gluons were tuned to OPAL LEP 1 data as described in [25,26]. We used 
HERWIG version 5.9 [20], PHOJET version 1.05c [27,28] and VERMASEREN version 1.01 [29] for 
two-photon interactions and KORALZ version 4.02 [30] for e^e^ — > t'^t" events. Generated 
events were processed through a full simulation of the OPAL detector [31] and the same 
event analysis chain was applied to the simulated events as to the data. 4,000,000 fully 
simulated events were generated by JETSET, 200,000 events, 1,000,000 events and 55,000 
events were generated by HERWIG, PHOJET and VERMASEREN while 800,000 events were 
generated by KORALZ. 

4 Event Selection 

4.1 Hadronic Event Selection 

This study is based on a sample of 3 million hadronic Z° decays selected as described 
in [32] from the data accumulated between 1992 and 1995 at cms energy of 91.2 GeV. 
We required that the central detector and the electromagnetic calorimeter were fully 
operational. 

For this study, we apply stringent cuts on tracks and clusters and further cuts on 
hadronic events. The clusters in the electromagnetic calorimeter are required to have 
a minimum energy of 100 MeV in the barrel and 250 MeV in the endcap. Tracks are 



required to have transverse momentum pt > 150 MeV/c with respect to the beam axis, 
at least 40 reconstructed points in the jet chamber, at the point of closest approach a 
distance between the track and the nominal vertex do < 2 cm in the r-0 plane and 
Zq < 25 cm in the z direction. We require at least five such tracks to reduce background 
from e+e" t^t~ and 77 qq events. The polar angle of the thrust axis is required 
to satisfy | cos6't| < 0.9, to ensure that events are well contained in the OPAL detector. 
After these cuts, a data sample of 2.4 x 10^ events remains. 



4.2 Isolated Photon Selection 
4.2.1 Isolation Cuts 

Isolated photons are selected in these hadronic events as follows. Electromagnetic clusters 
with an energy E-^c > 10 GeV are chosen in order to suppress background from soft 
photons coming from the decay of mesons. Accordingly, our signal event is defined as 
an e''"e~ — > qq event with an ISR or FSR photon with energy greater than 10 GeV. We 
use electromagnetic clusters in the polar angle region |cos^ec| < 0-72 corresponding to 
the barrel of the detector, where there is the least material in front of the lead glass, see 
Figure 2 a). Also, the non-pointing geometry of the endcap electromagnetic calorimeter 
complicates the cluster shape fitting explained below. The number of clusters in the data 
which satisfy Eec > 10 GeV and | cos6Iec| < 0.72 is 1,797,532. According to the Monte 
Carlo simulation, 99.3% of these selected clusters come from non-radiative multi-hadronic 
events. 

The candidate clusters are required to be isolated from any jets, and from other clusters 
and tracks: 

• The angle with respect to the axis of any jet, ctiso, is required to be larger than 25°, 
see Figure 2 b). The jets are reconstructed from tracks and electromagnetic clusters, 
excluding the candidate cluster, using the Durham algorithm [33] with ^cut = 0.005. 

• The sum of the momenta Piso of tracks falling on the calorimeter surface inside a 0.2 
radian cone around the photon candidate is required to be smaller than 0.5 GeV/c 
(Figure 2c)). The total energy deposition in the electromagnetic calorimeter within 
a cone of 0.2 radian around the photon candidate, £'iso, is also required to be less 
than 0.5 GeV (Figure 2 d)). 

After the isolation cuts, 11,265 clusters are retained. The fraction of clusters from non- 
radiative multi-hadronic events is reduced to 52.8%. The background from t~^t~ events 
(two-photon events) is 0.6% (0.01%) [34]. 



4.2.2 Likelihood Photon Selection 

Isolated photon candidates are selected by using a likelihood ratio method with four input 
variables, see appendix A for details. The first two variables are | cos^ecI and ojiso, defined 
above. Two more variables, the cluster shape fit variable S and the distance A between 
the electromagnetic calorimeter cluster and the associated presampler cluster, defined as 
follows, reduce the background from clusters arising from the decays of neutral hadrons. 
The cluster shape fit variable, S, is defined by 

Q 1 y^ (-E'meaSjj -^exp,i) ^g>^ 



-^block j "^measji 



where A^biock is the number of lead glass blocks included in the electromagnetic cluster, 
Emeas,i is the measured energy deposit in the ith block, -Eexp,i is the expected energy 
deposit in the ith block, assuming that the energy is deposited by a single photon, and 
c"meas,i is the Uncertainty in the energy measured by ith electromagnetic calorimeter block. 
Ecxp,i is a function of position and energy of the incident photon based on the simulation 
of the OPAL detector with single photons. The value of S is determined by minimizing 
Equation (6) under variation of the position and energy of the cluster. For a cluster to be 
considered further in the likelihood, preselection cuts are applied: we require the number 
of blocks to be at least two and the value of S after the fit to be smaller than 10. The 
quality of the cluster shape fits depends on the assumed resolution (Jmeas/, this will be 
studied as a systematic uncertainty. 

The variable A measures the distance between the electromagnetic calorimeter cluster 
and the associated presampler cluster, A = max(|A0|, |A^|), with A0 and A^ the angular 
separations between the clusters. 

The distributions of 5* and A are shown in Figures 2 e) and f). The Monte Carlo 
distributions in these figures are normalized according to the luminosity obtained from 
small angle Bhabha events. 

A disagreement between data and Monte Carlo is seen for S and ctiso. The level of 
agreement between data and Monte Carlo for the S distribution is studied with photons 
in radiative muon pair events and tt^s produced in tau pair events. It is confirmed that the 
Monte Carlo adequately reproduces the S distributions [34]. The disagreement between 
data and Monte Carlo for distributions of S and aiso stems from the failure of the Monte 
Carlo generators to correctly predict the rate of isolated neutral hadrons, as explained in 
Section 4.4. In this analysis, the rate of isolated neutral hadrons used in the background 
subtraction is estimated from data by methods described in Section 4.4. 

The likelihood calculation is performed with reference histograms made for seven sub- 
samples, chosen according to the cluster energy. The cut on the likelihood value is chosen 
so as to retain 80% of the signal events. The hkelihood distributions for data and Monte 
Carlo arc shown in Figure 3. It can be seen that the likelihood distributions for signal and 
background events are well separated for each region of electromagnetic cluster energy. 
Electromagnetic clusters which pass the likelihood selection are regarded as photon can- 
didates. If more than one candidate is found in the same event the one with the highest 
energy is chosen. 

4.3 Final Data Sample 

Hadronic events with hard isolated photon candidates are divided into seven subsamples 
according to the photon energy for further analysis. Table 1 shows the mean values of 
Vs', the number of data events and the number of background events for each subsample. 

4.4 Background Estimation 

According to the Monte Carlo simulation, the contamination from r pair events is between 
0.5 and 1.0%. The impact of this small number of events is further reduced because the 
value of event shape observables for r pair events are concentrated in the lowest bin of 
the distributions, outside the fitting range, so their effect on the fits is negligible. The 
contribution of two photon processes is less than 0.01% in all subsamples and is ignored. 



As mentioned in [35], the JETSET Monte Carlo fails to reproduce the observed rate 
of isolated electromagnetic clusters, both for isolated photons and isolated 7r°'s. Isolated 
neutral hadrons are the dominant source of background for this analysis, and their rate 
has been estimated from data using the following two methods. 

Firstly, with the likelihood ratio method the observed likelihood distributions in the 
data in bins of photon energy were fitted with a linear combination of the Monte Carlo 
distributions for signal and background events which pass the isolation cuts and hkelihood 
preselection requirements. The overall normalisation of the Monte Carlo distribution is 
fixed to the number of data events. The fit uses a binned maximum likelihood method 
with only the fraction of background events as a free parameter. Figure 3 shows the fit 
results. The values of x^/d.o.f. are between 1.2 and 3.4 for 18 degrees of freedom. 

Secondly, with the isolated tracks method the fraction of background from isolated 
neutral hadrons was estimated from the rates of isolated charged hadrons. We select from 
the data tracks which satisfy the same isolation criteria as the photon candidates. The 
composition of these isolated charged hadrons obtained from JETSET is used to infer the 
rates of charged pions, kaons and protons. When isospin symmetry is assumed, the rates 
of neutral pions, neutral kaons and neutrons can be estimated from the rates of charged 
pions, charged kaons and protons, respectively: 

-RttO — o-^Tri Rk° — o-Rk± Rn — Rp , (7) 

z ^ z 

where Rx is the production rate of particle X. According to JETSET tuned with OPAL 
data, the rate of isospin symmetry violation is 10% for pions and 5% for kaons and protons. 
This is assigned as a systematic uncertainty for the isolated tracks method and combined 
with the statistical uncertainty. 

The neutral hadron background fractions estimated by these two methods are shown 
in Table 1. The statistical errors from the number of data and Monte Carlo events from 
fitting the hkelihood distributions are shown. The results from the two methods are within 
at most three standard deviations of these errors, except in the bin 35 — 40 GeV. 

The standard analysis will use the likelihood ratio method. Any differences in the 
resulting values of Og obtained by using the two background estimate methods will be 
treated as a systematic uncertainty. 

5 Measurement of Event Shape Distributions 

In this analysis event shape observables as defined above in section 2 are calculated 
from tracks and electromagnetic clusters excluding the isolated photon candidate. The 
contributions of electromagnetic clusters originating from charged particles are removed 
by the method described in [36]. 

We evaluate the observables in the cms frame of the hadronic system. The Lorentz 
boost is determined from the energy and angle of the photon candidate. When the four- 
momentum of particles in the hadronic system is calculated, electromagnetic clusters are 
treated as photons with zero mass while tracks of charged particles are treated as hadrons 
with the charged pion mass. 

Distributions of the event shape observables (1 — T) and Mh are shown for two cms 
energies in Figure 4. The remaining background is removed by subtracting the scaled 
Monte Carlo predictions for non-radiative hadronic events and for r pair events using the 



background estimates listed in Table 1. The effects of the experimental resolution and 
acceptance are unfolded using Monte Carlo samples with full detector simulation (detector 
correction). The unfolding is performed bin- by-bin with correction factors rf^* = hi/di, 
where hi represents the value in the ith bin of the event shape distribution of stable 
hadrons in the Monte Carlo simulation, where "hadrons" are defined as particles with a 
mean proper lifetime longer than 3 • 10~^°s. di represents the value in the ith bin of the 
event shape distribution calculated with clusters and tracks obtained from Monte Carlo 
samples with detector simulation after the complete event selection has been applied. We 
refer to the distributions after applying these corrections as data corrected to the hadron 
level. 

The distributions of the event shape observables 1 — T and Mh for data corrected to 
the hadron level and corresponding Monte Carlo predictions are shown in Figure 5^. The 
Monte Carlo samples are generated with cms energies set to the mean value of in 
each subsample. In the production of the Monte Carlo samples ISR and FSR is switched 
off and on, respectively. The predictions from the event generators are consistent with 
the data for all x/i' bins. There is similar agreement between data and event generator 
predictions for the other observables. 

6 Measurement of cxs 

The measurement of is performed by fitting perturbative QCD predictions to the event 
shape distributions corrected to the hadron level for (1 — T), Mh [37], Brp^ By/ [38,39], 
C [40,41] and yg [33,42,43]. The 0{a^) and NLLA calculations are combined with the 
ln(i?) matching scheme. The effects of hadronisation on event shapes must be taken into 
account in order to perform fitting at the hadron level (hadronisation correction). Pre- 
serving the normalisation in the hadronisation correction is not trivial for low x/s' samples 
because of large hadronisation corrections. The hadronisation correction is applied to the 
integrated (cumulative) theoretical calculation to conserve normalisation as in our pre- 
vious analyses [44,45]. The hadron level predictions are obtained from the cumulative 
theoretical calculation multiplied by a correction factor Rf"^ — Hi/ Pi, where Pi (Hi) 
represents the value in the ith bin of the cumulative event shape distribution calculated 
by Monte Carlo simulation without (with) hadronisation. The JETSET Monte Carlo event 
generator is used for our central results, while HERWIG and ARIADNE are considered as 
alternatives for the estimation of systematic uncertainties. 

The fit of the hadron level QCD predictions to the event shape observables uses a 
least method with as{Q) treated as a free parameter. Only statistical uncertainties 
are taken into account in the calculation of x^. When the total number of events is small, 
the differences between the statistical errors counting larger or smaller numbers of events 
than the theoretical prediction can bias the fit result. In order to avoid this bias the value 
of the fitted theoretical distribution is used to calculate the statistical error instead of 
the number of events in each bin of the data distribution. The statistical uncertainty is 
estimated from the fit results derived from 100 Monte Carlo subsamples with the same 
number of events as selected data events. 

The region used in the fit is adjusted such that the background subtraction and the 
detector and hadronisation corrections are small (less than 50%) and uniform in that 

^The values of the six observables at the seven energy points are given in [34] and will be available 
under http : //durpdg . dur . ac . uk/HEPDATA/. 



region. The resulting fit ranges are mainly restricted by the hadronisation corrections. 
The QCD predictions at y/s' — 78 GeV fitted to data after applying the hadronisation 
correction are shown in Figure 6. Good agreement between data and theory is seen. The 
fitted values of ag and their errors for each event shape observable are shown in Tables 2-5. 

6.1 Systematic Uncertainties 

6.1.1 Experimental Uncertainties 

The experimental uncertainty is estimated by adding in quadrature the following contri- 
butions: 

• the difference between the standard result and the result when all clusters and tracks 
are used without correcting for double counting of energy. This variation is sensitive 
to imperfections of the detector simulation. 

• the largest deviation between the standard result and the result when the analysis is 
repeated with tighter selection criteria to eliminate background (standard values in 
brackets): the thrust axis is required to lie in the range | cos^t| < 0.7 (0.9), or the 
cluster shape variable is required to be smaller than 5 (10), or the isolation angle 
from any jet is required to be larger than 35° (25°). 

• the difference between the standard value and the value obtained by repeating the 
analysis with the background fractions estimated from the rate of isolated charged 
hadrons as described in Section 4. 

• the difference between the standard result and the result when the single block 
energy resolution is varied to give the expected the cluster shape fits. This 
check is made, because the values of in the cluster shape fits depend on the 
assumed energy resolution. 

• the maximum difference between the standard result and the result when the fit 
regions are varied. The lower and upper limit of the fitting region are independently 
changed by ±1 bin. 

The tighter selection on | cos^r| and the alternative single block energy resolution of the 
electromagnetic calorimeter yield the largest contributions to the experimental systematic 
uncertainty. The overall resolution and energy scale uncertainty of the electromagnetic 
calorimeter have a neglegible effect on the results of this analysis. 

6.1.2 Hadronisation Uncertainties 

The following variations are performed in order to estimate the hadronisation uncertain- 
ties: 

• the largest of the changes in ag observed when independently varying the hadroni- 
sation parameters h and ctq by ±1 standard deviation about their tuned values in 
JETSET [25]; 

• the change observed when the parton virtuality cut-off parameter Qq is varied by 
±1 standard deviation about its tuned value in JETSET; 



• the change observed when only the hght quarks u, d, s and c are considered at the 
parton level in order to estimate potential quark mass effects; 

• the differences with respect to the standard result when HERWIG or ARIADNE are used 
for the hadronisation correction, rather than JETSET. 

We define the hadronisation correction uncertainty by adding in quadrature the deviation 
when using only light quarks and the larger deviation when using HERWIG or ARIADNE to 
calculate the corrections. These variations are observed to lead to larger differences than 
all other variations, i.e. the main contributions to the hadronisation uncertainties are the 
choice of hadronisation model and the potential effect of quark masses. 

6.1.3 Theoretical Uncertainties 

We fix the renormalisation scale parameter = yu/Q to 1, where fi is the energy scale at 
which the theory is renormalized and Q is the energy scale of the reaction. Although the 
uncertainty on the choice of the value of X/^ gives a large contribution to the systematic 
uncertainty, the means of quantifying this uncertainty is essentially arbitrary. We define 
the scale uncertainty as the larger of the deviations of as when is changed from 1 to 
0.5 or 2.0. 

The 0{al) and NLLA calculations are combined with the ln(R) matching scheme. 
The variation in ^^(v^) due to using different matching schemes is much smaller than 
the renormalisation scale uncertainty [46] , and is not included as an additional theoretical 
systematic uncertainty. 

6.2 Combination of Results 

The values of ag obtained for each observable at each energy are used to study the energy 
dependence of ag and to obtain an overall combined result for as{Mz)- The individual 
values of a.^ as given in Tables 2-5 and shown in Figure 7 are combined taking the cor- 
relations between their statistical and systematic errors into account using the method 
described in [8] . The statistical covariances between results from different observables are 
determined at each energy from 100 Monte Carlo subsamples with the same number of 
events as selected in the data. The experimental systematic uncertainties are assumed 
to be partially correlated, i.e. covij = min((Tj, 0"^)^. The hadronisation and theoretical 
covariances are only added to the diagonal of the total covariance matrix. The correla- 
tions between these uncertainties are considered by repeating the combination procedure 
with different hadronisation corrections (udsc only, HERWIG, ARIADNE) and with different 
renormalisation scale parameters {x^ = 2 and 0.5). The systematic uncertainties for the 
combined value are obtained by repeating the combination for each systematic variation. 
The resulting values of as{\/s') are shown in Table 6 and Figure 8. 

Values of as from individual observables at each energy are combined after evolving 
them to a/s — Mz- In this case the results are statistically uncorrelated. The correlations 
between systematic uncertainties are treated as explained above. The results are given in 
Table 7 and Figure 9. 

We also combine the combined values listed in Table 7 taking into account their 
statistical correlations using the sum of the inverses of the individual statistical covariance 
matrices at each energy point. The result is 



as{Mz) = 0.1182 ± 0.0015(stat.) ± 0.0101(syst.). 



(8) 



and is shown with individual errors in Table 7. Figure 8 shows the evolution of the 
strong coupling using our result. As a crosscheck on the robustness of the combination 
procedure we repeat the combination using the combined results at each energy point 
shown in Tabic 6 or using all individual results and find as{Mz) — 0.1183 ± 0.0103 or 
a^(Mz) = 0.1179 ± 0.0103, respectively. 

Our result is consistent within the statistical and experimental errors with the result 
from OPAL using non-radiative events in LEP 1 data with the same set of observables, 
as{Mz) = 0.1192±0.0002(stat.)±0.0050(syst.) [8]. Our result is also consistent with recent 
combined values [11,47-49] and results from other analyses using radiative events [16,17]. 
Figure 10 compares our result with the LEP 1 value for as from [8] and an average of 
results from L3 using radiative hadronic Z° decays [16]^. 

The combinations of individual observables at different cms energies yield x^/d.o.f. ^ 
1/6. The small values of x^/d.o.f. are due to the conservative treatment of hadronisation 
and theoretical uncertainties. The values of x^/d.o.f. indicate consistency of the individual 
results with the model of the combination including evolution of results at different cms 
energies to Mz before the combination. 

7 Summary 

The strong couphng ag has been measured at reduced cms energies, y/s', ranging from 
20 GeV to 80 GeV using event shape observables derived from the hadronic system in 
radiative hadronic events. 

Fits of 0{af) and NLLA QCD predictions to the six event shape observables 1 — T, 
Mh, Bt, BiY, C and |/£ are performed and values of are obtained for seven values of 
sfs'. Our results are consistent with the running of ag as predicted by QCD. The values 
at each y/s' are evolved to jj, — Mz and combined for each event shape observable. The 
combined value from all event shape observables and values is as{Mz) — 0.1182 ± 
0.0015(stat.) ±0.0101(syst.). 

This result agrees with previous OPAL analyses with non-radiative LEP 1 data, with 
a similar measurement by L3, and with recent world average values, see figure 10. Within 
errors, QCD is consistent with our data sample of events with isolated FSR. Our result 
supports the assumption that the effects of high energy and large angle FSR on event 
shapes in hadronic Z° decays can be effectively described by QCD with a lower effective 
cms energy \fs' . 
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A Likelihood Ratio Method 

The likelihood ratio C^q^ is defined by 

r - = (Q) 

where L^q^ and Lbkg,i are the absolute likelihood values for signal qq7 events and events 
from the ith background process. The background likelihood values are weighted by Wi 
proportional to the cross section of the ith background process. 

The absolute likehood values L are calculated from probability density functions (pdfs) 
Pj{xj) for the input variables Xj. The pdfs arc obtained as so-called reference histograms 
from simulated signal and background samples. For the calculation of the pdfs the pro- 
jection and correlation approximation (PCA) method [50] is used. In brief, each Xj is 
transformed to a variable yj following a Gaussian distribution using 

y = y/2eii-\2F{x) - 1) (10) 

where erf~^ is the inverse error function and F{x) — J^^.^p{x')dx' is the cumulative 
distribution of x. The hkelihood L{x) is then given in the PCA by 

m = -^e-y^^'^-'-'^y/'Up^ix,) . (11) 



V is the n X n covariance matrix of the yj, I is the identity matrix and x and y are the 
vectors of the Xj and yj. 
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Table 1: The number of selected events and the mean value of y/s' for each y/s' subsample. 
The neutral hadron background fractions estimated by the two methods described in 
Section 4.4 are listed in the columns "Non-rad. MH". 
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Table 2: Values of as and their errors for subsamples £'.^=10-15 GeV (upper) and 15- 
20 GeV (lower). 
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-0.0064 


-0.0014 


x^, = 2.0 


0.0073 


0.0058 


0.0076 


0.0048 


0.0081 


0.0048 


Total error 


+ 0.0136 
- 0.0128 


+ 0.0117 
- 0.0111 


+ 0.0134 
- 0.0126 


+ 0.0102 
- 0.0096 


+ 0.0148 
- 0.0140 


+ 0.0195 
- 0.0190 


(1 - T) Mh 


Bx 


Bw 


C 




as(57.6GeV) 0.1378 


0.1396 


0.1327 


0.1194 


0.1284 


0.1407 


Sta.tisticiil Error 


± O.OOS-l 


± 0.0094 


+ (1 (1(179 








Tracks + Clusters 


0.0004 


0.0022 


—0.0008 


0.0005 


n nn o n 

0.0039 


—0.0013 


|cos6't| < 0.7 


0.0065 


0.0101 


0.0078 


0.0054 


0.0083 


0.0056 


C>5 


-0.0003 


0.0020 


U.UUio 


n nni o 
U.UUiC) 


n nnnc 
U.UUUo 


n nnno 

U.uuuy 


af° 


-0.0010 


-0.0052 


0.0004 


-0.0004 


0.0001 


-0.0007 


Bkg fraction 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


ECAL Resolution 


0.0032 


0.0051 


0.0035 


0.0021 


0.0010 


-0.0013 


Fitting Range 


0.0036 


0.0006 


0.0014 


0.0020 


0.0011 


0.0010 


Experimental Syst. 


± 0.0082 


± 0.0128 


-1- n nnss 


m u.uuoo 


m u.uuyo 


m u.uuoi 


6- Is.d. 


-0.0009 


-0.0004 


-0.0005 


-0.0003 


-0.0010 


-0.0006 


6+ Is.d. 


0.0006 


0.0004 


0.0004 


0.0004 


0.0009 


0.0005 


Qo - Is.d. 


0.0002 


-0.0008 


0.0006 


-0.0005 


0.0005 


-0.0023 


Qo + Is.d. 


-0.0005 


0.0005 


—0.0009 


0.0003 


—0.0004 


0.0016 


aq — Is.d. 


0.0005 


0.0002 


O.UOOo 


U.UUUo 


u.uuuy 


U.UUUo 


aq + Is.d. 


-0.0009 


-0.0002 


—0 0005 


—0.0004 


—0.0010 


—0.0007 


udsc only 


0.0024 


-0.0001 


0.0042 


0.0033 


0.0040 


0.0063 


Iforwig 5.9 


-0.0076 


-0.0039 


-0.0072 


-0.0066 


-0.0101 


-0.0113 


Ariadne 4.08 


-0.0011 


-0.0011 


-0.0012 


-0.0013 


-0.0032 


-0.0049 


Total Hadronisation 


± 0.0081 


± 0.0041 


± 0.0085 


± 0.0075 


± 0.0114 


± 0.0140 


a;^ = 0.5 


-0.0079 


-0.0063 


-0.0078 


-0.0042 


-0.0072 


-0.0023 


a;^ = 2.0 


0.0101 


0.0087 


0.0098 


0.0058 


0.0090 


0.0063 


Total error 


+ 0.0175 
- 0.0164 


+ 0.0186 
- 0.0176 


+ 0.0172 
- 0.0162 


+ 0.0130 
- 0.0124 


+ 0.0183 
- 0.0175 


+ 0.0189 
- 0.0180 



Table 3: Values of as and their errors for subsamples — 20 — 25 GeV (upper) and 
25 -30 GeV (lower). 



(1 - T) Mh 




Bw 


C 




Q!^(49.0GeV) 0.1373 0.1359 


0.1413 


0.1269 


0.1356 


0.1440 


Statistical Error 


± 0.0105 


± 0.0098 


HZ U.UUO ( 


m u.uuuy 


-1- n nnsQ 


-1- n 01 1 7 


Tracks + Clusters 


0.0022 


0.0007 


0.0032 


r\ Anno 

—0.0003 


n nnn o 

0.0008 


—0.0012 


cos6't < 0.7 


0.0029 


0.0039 


r» AAA A 

0.0004 


0.0012 


—0.0001 


A AAAA 

—0.0000 


C > 5 


—0.0010 


—0.0038 


— U.UUi ( 


— U.UUUi 


— u.uu4y 


n nm q 

— U.UUio 




0.0024 


0.0024 


0.0007 


0.0017 


0.0013 


0.0046 


Bkg fraction 


0.0001 


0.0000 


0.0001 


0.0000 


0.0000 


0.0001 


ECAL Resolution 


-0.0003 


0.0010 


-0.0003 


0.0009 


-0.0000 


-0.0005 


Fitting Range 


0.0027 


0.0013 


0.0009 


0.0016 


0.0009 


0.0020 


Experimental Syst. 


± 0.0053 


± 0.0062 


m u.uuoo 


-1- n nn9s 
m u.uuzo 


m u.uuoz 


m u.uuoo 


h-ls.d. 


-0.0005 


-0.0009 


—0.0006 


—0.0005 


—0.0009 


—0.0008 


b+ Is.d. 


0.0005 


0.0008 


0.0003 


0.0005 


0.0007 


0.0002 


Qa - Is.d. 


0.0003 


-0.0006 


0.0006 


-0.0003 


0.0005 


-0.0019 


Qo + Is.d. 


-0.0005 


0.0005 


-0.0012 


0.0004 


-0.0005 


0.0017 


aq — Is.d. 


0.0005 


0.0005 


U.UUU4 


U.UUU 1 


r\ nnn'7 
U.UUU ( 


n nnnT 
U.UUU/ 


(Tq + is.d. 


— U.OUUd 


— 0.0006 


—0 0008 


—0 0007 


— OOOQ 

yj .\J\J\JiJ 


— OOOfi 


udsc only 


0.0023 


0.0002 


0.0039 


0.0050 


0.0038 


0.0060 


TTpfixri cr ^ Q 


_n 0083 


\J •\J\JiJ\J 


-0.0080 


-0.0083 


-0.0123 


-0.0114 


Ariadne 4.08 


-0.0009 


-0.0041 


-0.0011 


-0.0024 


-0.0039 


-0.0056 


Total Hadronisation 


± 0.0087 


± 0.0099 


± 0.0091 


± 0.0101 


± 0.0135 


± 0.0142 


= 0.5 


-0.0076 


-0.0058 


-0.0092 


-0.0054 


-0.0081 


-0.0008 


x^, = 2.0 


0.0097 


0.0081 


0.0117 


0.0072 


0.0102 


0.0056 


Total error 


+ 0.0176 
- 0.0165 


+ 0.0173 
- 0.0163 


+ 0.0176 
- 0.0160 


+ 0.0144 
- 0.0136 


+ 0.0198 
- 0.0188 


+ 0.0201 
- 0.0193 


(1 - T) Mh 


Bx 


Bw 


C 




as(38.5GeV) 0.1474 


0.1374 


0.1451 


0.1415 


0.1421 


0.1496 


Sta.tisticiil Error 


± O.OVir, 


± 0.0112 


□Z \> .KiKjOO 


+ (1(111'-! 


+ (1 oil'-! 




Tracks + Clusters 


0.0024 


0.0019 


0.0006 


0.0001 


0.0049 


—0.0010 


|cos6't| < 0.7 


0.0026 


0.0059 


0.0034 


0.0061 


0.0050 


A AAAA 

0.0022 


C>5 


0.0042 


0.0038 


U.UUio 


U.UUo ( 


U.UUOZ 


U.UU4U 


af° 


0.0005 


-0.0007 


-0.0004 


0.0043 


0.0014 


0026 


Bkg fraction 


0.0003 


0.0003 


0.0002 


0.0002 


0.0003 


0.0004 


ECAL Resolution 


0.0019 


0.0025 


0.0003 


0.0035 


0.0039 


0.0055 


Fitting Range 


0.0033 


0.0009 


0.0008 


0.0013 


0.0023 


0.0008 


Experimental Syst. 


± 0.0067 


± 0.0077 








-1- n 0077 
m u.uu ( ( 


6- Is.d. 


-0.0009 


-0.0007 


-0.0007 


-0.0004 


-0.0007 


-0.0005 


6+ Is.d. 


0.0009 


0.0006 


0.0006 


0.0005 


0.0005 


0.0004 


Qo - Is.d. 


0.0006 


-0.0008 


0.0011 


-0.0008 


0.0007 


-0.0021 


Qo + Is.d. 


-0.0005 


0.0008 


—0.0014 


0.0006 


—0.0007 


0.0018 


aq — Is.d. 


0.0013 


0.0003 


O.UOOo 


U.OUUd 


n nni n 
U.UUiU 


n nnnc; 
O.UUUO 


aq + Is.d. 


-0.0009 


-0.0002 


—0 0008 


—0.0004 


—0.0007 


—0 OOOfi 


udsc only 


0.0042 


0.0001 


0.0060 


0.0036 


0.0038 


0.0064 


Iforwig 5.9 


-0.0150 


-0.0096 


-0.0105 


-0.0107 


-0.0125 


-0.0127 


Ariadne 4.08 


-0.0042 


-0.0036 


-0.0028 


-0.0025 


-0.0030 


-0.0055 


Total Hadronisation 


± 0.0162 


± 0.0103 


± 0.0125 


± 0.0116 


± 0.0135 


± 0.0154 


a;^ = 0.5 


-0.0093 


-0.0055 


-0.0097 


-0.0072 


-0.0089 


-0.0012 


a;^ = 2.0 


0.0120 


0.0079 


0.0124 


0.0097 


0.0114 


0.0063 


Total error 


+ 0.0247 
- 0.0235 


+ 0.0188 
- 0.0179 


+ 0.0201 
- 0.0186 


+ 0.0210 
- 0.0199 


+ 0.0232 
- 0.0221 


+ 0.0210 
- 0.0200 



Table 4: Values of as and their errors for subsamples — 30 — 35 GeV (upper) and 
35 -40 GeV (lower). 



(1 - T) Mh 






C 




Q!^(24.4GeV) 0.1569 0.1524 


n 1 

u. ±oo^ 


U. ±'-tOO 


14.0fi 




Statistical Error 


± 0.0252 


± 0.0117 


± 0.0115 


± 0.0101 


± 0.0112 


± 0.0181 


Tracks + Clusters 


0.0038 


0.0015 


0.0060 


-0.0021 


0.0080 


-0.0074 


Icos(?rI < 0.7 


0.0001 


0.0008 


-0.0027 


0.0027 


0.0013 


-0.0008 


C > 5 


0.0037 


-O.OOOl 


-0.0U3G 


-0.0022 


-0.0010 


-0.0084 


ar 


0.0110 


0.0056 


0.0003 


0.0023 


0.0060 


0.0005 


Bkg fraction 


0.0023 


0.0017 


n An 1 o 

0.0018 


A AA1 C 

0.0015 


A A AO A 

0.0020 


A A AO 1 

0.0031 


ECAL Resolution 


-0.0035 


-0.0053 


—0.0039 


—0.0013 


—0.0025 


—0.0057 


Fitting Range 


0.0035 


0.0027 


0.0018 


A AA1 'T 

0.0017 


A A AO A 

0.0020 


A A A 1 O 

0.0018 


Experimental Syst. 


± 0.0134 


± 0.0085 


± 0.0088 


± 0.0054 


± 0.0109 


± 0.0131 


b-ls.d. 


-0.0007 


-0.0014 


-0.0007 


-0.0013 


-0.0012 


-0.0013 


b+ls.d. 


0.0015 


0.0017 


0.0009 


0.0011 


0.0012 


0.0006 


Qo - Is.d. 


0.0010 


-0.0010 


0.0023 


-0.0009 


0.0010 


-0.0039 


Qo + Is.d. 


-0.0008 


0.0004 


-0.0029 


0.0000 


-0.0010 


0.0015 


(Jq — Is.d. 


0.0014 


0.0010 


0.0011 


0.0017 


0.0012 


0.0011 


Uq + Is.d. 


-0.0010 


-0.0009 


-0.0010 


-0.0018 


-0.0013 


-0.0016 


udsc only 


0.0075 


0.0053 


0.0140 


0.0159 


0.0150 


0.0168 


TT • r f\ 

Her wig 5.9 


—0.0212 


—0.0080 


— U.Uio4 


— U.UiZD 


— U.UiUo 


n ni QQ 
— u.uiyo 


Ariadne 4.08 


-0.0082 


-0.0056 


-0.0040 


-0.0045 


-0.0050 


-0.0114 


Total Hadronisation 


± 0.0240 


± 0.0113 


-1- n 0900 


-1- 0900 


+ 01 QO 
□z u.uiyu 


-1- 098'? 


Xn = 0.5 


-0.0104 


-0.0085 


-0.0116 


-0.0082 


-0.0088 


-0.0024 


x^ = 2.0 


0.0137 


0.0115 


0.0151 


0.0108 


0.0112 


0.0084 


Total error 


+ 0.0397 
- 0.0387 


+ 0.0216 
- 0.0202 


+ 0.0289 
- 0.0273 


+ 0.0262 
- 0.0252 


+ 0.0270 
- 0.0261 


+ 0.0371 
- 0.0362 



Table 5: Values of ag and their errors for subsample — AO — 45 GeV. 



Vs' [GeV] 


78.1 71.8 65.1 57.6 49.0 38.5 24.4 




0.1153 0.1242 0.1201 0.1296 0.1353 0.1438 0.1496 


Statistical 
Experimental 
Hadronisation 
Theory 


0.0026 0.0037 0.0039 0.0047 0.0053 0.0064 0.0071 
0.0068 0.0036 0.0040 0.0069 0.0039 0.0063 0.0077 
0.0062 0.0065 0.0072 0.0085 0.0100 0.0122 0.0166 
0.0053 0.0067 0.0063 0.0076 0.0086 0.0099 0.0117 



Table 6: Combined values of as{\/s') and their errors from all event shape variables. 





(1-T) Mh Bt Bw C y^^ 


Combined 


asiMz) 


0.1230 0.1187 0.1214 0.1117 0.1195 0.1261 


0.1182 


Statistical 
Experimental 
Hadronisation 
Theory 


0.0028 0.0024 0.0021 0.0021 0.0023 0.0031 
0.0050 0.0054 0.0037 0.0033 0.0040 0.0049 
0.0071 0.0052 0.0080 0.0061 0.0092 0.0105 
0.0076 0.0059 0.0081 0.0049 0.0076 0.0045 


0.0015 
0.0038 
0.0070 
0.0062 



Table 7: Combined values of as{Mi) and their errors from all photon energy subsamples 
for a given observable. The final combined value of as{Mz) is also shown. 






Figure 1: The distributions of event shape observables 1 — T and Mh for non-radiative 
events and radiative hadronic events from the Monte Carfo generators JETSET, HERWIG 
and ARIADNE as indicated befow the figures. The triangles and points show distributions 
obtained from the Z° samples with FSR while the histograms show distributions from 
samples generated at lower energies as shown on the figure. The open triangles and solid 
histogram (solid points and dashed histogram) in each figure correspond to y/s' = 40 
(70) GeV. 
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Figure 2: Distributions of each variable used in the isolated photon selection. The error 
bars show the statistical errors. Monte Carlo distributions are normalized to the integrated 
luminosity of the data and the cross section of the process. Arrows in the figures show the 
selected region. Distributions for radiative multi-hadronic events, which are signal events 
in this analysis, are overlaid on distributions for all multi-hadronic events and rr events. 
The distribution of each variable is obtained with the cuts on the preceeding variables 
applied. 
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Figure 3: Photon likelihood distributions. The error bar shows statistical error. The 
Monte Carlo distributions are normalised to the total number of candidates in the data, 
and the neutral hadron background fractions are obtained from the fits described in 
Section 4.4. The arrows indicate the selected regions. 
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Figure 4: Event shape distributions before background subtraction and detector correc- 
tion. Two of the six event shape observables, 1 — T and Mh, are shown for the low 
(38.5 GeV) and high (78.1 GeV) Vs' samples. The histograms show Monte Carlo distri- 
butions. The error bars show the statistical errors. 
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Figure 5: Event shape distributions at the hadron level. The error bars correspond to 
the statistical and experimental uncertainties described in Section 6.1.1. Two of the 
six event shape observables, 1 — T and M/^, are shown for the low (38.5 GeV) and high 
(78.1 GeV) \fs' samples. The small lines on the error bars show the extent of the statistical 
uncertainty. The data points are placed at the centres of the corresponding bins. The 
predictions of JETSET, HERWIG and ARIADNE at the corresponding values are also 
shown as lines. 
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Figure 6: Event shape distributions for data at (v s') = 78.1 GeV and the fitted theoretical 
predictions. The error bars show the statistical errors. The solid lines in the theoretical 
predictions show the regions used in the fit. Three corrections are plotted as "i?corr": 
the detector correction, rf^* (dashed line), the hadronisation correction, Rf"''^ (solid line), 
and the ratio of distributions after and before background subtraction (dotted line). The 
hadronisation correction is shown by the ratio of differential distributions in these figures 
(see text for details). 
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Figure 7: Energy dependence of as for all ^fs' subsamples. The inner error bars show the 
statistical and the outer error bars the total uncertainties. The curves and shaded bands 
show the QCD prediction for the running of as obtained with the corresponding values 
of as{Mi) with total errors from Table 7. 




Figure 8: Combined values of from all event shape observables as shown in Table 6. 
The curve and shaded band show the QCD prediction for the running of as using the 
combined value of as{Mz) with total errors. 
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Figure 9: The values of as{Mz) obtained by combining all y/s' samples as shown in 
Table 7. The inner error bars are statistical, the outer error bars correspond to the total 
uncertainty. The dashed vertical lines and shaded bands show the LEP 1 results from 
OPAL [8] using non-radiative events. 
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Figure 10: Combined values of as{Mz) for all event shape observables and \/^ samples. 
The error bars show total uncertainties. The results from this analyses with radiative 
events, from non-radiative events with OPAL LEP 1 data [8] and from L3 radiative 
events [16] are shown. The PDG [47] value of as{Mz) is also shown as the vertical line, 
with the total uncertainty corresponding to the shaded band. 



